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We investigate the role of electronic-vibrational coupling in resonant electron transport through 
single-molecule junctions, taking into account that the corresponding coupling strengths may de¬ 
pend on the charge and excitation state of the molecular bridge. In the presence of multiple elec¬ 
tronic states, this requires to extend the commonly used model and include vibrationally dependent 
electron-electron interaction. We use Born-Markov master equation methods and consider selected 
models to exemplify the effect of the additional interaction on the transport characteristics of a 
single-molecule junction. In particular, we show that it has a significant influence on local cool¬ 
ing and heating mechanisms, may result in negative differential resistance, and cause pronounced 
asymmetries in the conductance map of a single-molecule junction. 


I. INTRODUCTION 


Electron transport through a single-molecule junc¬ 
tion constitutes a complex many-body problem. Elec¬ 
tronic and vibrational degrees of freedom of a molec¬ 
ular conductor are often strongly correlated, in 
particular in nonequilibrium states at higher bias 
voltagesic—. The investigation of coupling between 
electronic and nuclear degrees of freedom in nanos¬ 
tructures under nonequilibriuin conditions has been 
of great interest recently and revealed a wealth of 
physical phenomena such as rectification^ii^^— , vibra¬ 
tionally induced decoherence^^^— negative differential 
resistanc o^^’^^i^°~ — , and signatures of multistability^^— . 

Most theoretical studies of vibrationally coupled elec¬ 
tron transport in molecular junctions have employed 
a simplified model of a molecule with linear coupling 
of the electronic degrees of freedom to vibrational 
modes described in the harmonic approximatio n^i^^i^^~ — 
More realistic models include charge-dependent vi¬ 
brational frequencie a^^d^~ — or anharmonic nuclear 
potential a^^d^~ — . The description of systems with 
multiple electronic states often requires to include the 
Coulomb interaction between electrons. The correspond¬ 
ing interaction strengths depend on the specific electronic 
configuration and on the nuclear geometry. This results 
in vibrationally dependent electron-electron interactions. 
The study of this effect and its manifestation in transport 
characteristics of molecular junctions is the main subject 
of this article. The importance of vibrationally depen¬ 
dent electron-electron interactions has already been real¬ 
ized, for example, in the context of the Hubbard model 
approach to electron-electron interaction in molecules^ 
as well as in studies of dissociation in colloidal quantum 
dot systems^. 

The outline of this article is as follows: The theoret¬ 
ical methodology is introduced in Sec. im including the 
model Hamiltonian of the molecular junctions, the mas¬ 


ter equation approach used to describe charge transport, 
and a discussion of the relevance of vibrationally depen¬ 
dent electron-electron interactions. In Sec. imi we an¬ 
alyze the effects and manifestation of vibrationally de¬ 
pendent electron-electron interactions in transport char¬ 
acteristics. To this end, we consider selected model sys¬ 
tems, including scenarios with symmetric and asymmet¬ 
ric molecule-lead coupling. 


II. THEORETICAL METHODOLOGY 

A. Model Hamiltonian 


The Hamiltonian of the molecular bridge in atomic 
units is given by 
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where i?nuc describes the nuclear and i7ei(R) the elec¬ 
tronic degrees of freedom of the molecule. The nuclear 
part of the Hamiltonian, iJnuc, includes the kinetic en¬ 
ergy of the nuclei and the Coulomb repulsion, I4uc(R)- 
Similarly, i7ei(R) includes the kinetic energy of the elec¬ 
trons, the Coulomb repulsion, and the Coulomb attrac¬ 
tion between the electrons and the nuclei. Thereby, the 
vector R summarizes the coordinates Ra of the nuclei 
and Yi denotes the coordinates of the Hh electron. The 
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charge and the mass of the nuclei are given by Za and 
Mq, respectively. 

We follow the scheme outlined to derive an 

approximate representation of the Hamiltonian in second 
quantization. Thereby we focus on the molecular elec¬ 
tronic states, thus working in a restricted subspace only 
incorporating a subset of all possible states. Accordingly, 
the result will be a restricted model Hamiltonian, con¬ 
taining effective interaction terms. Those effective inter¬ 
action terms not only account for the bare interaction, 
but also for the influence of the states beyond the re¬ 
stricted subspace under consideration. As a reference 
state we use the electronic ground state of the uncharged 
molecule 

77el(R) |«'ref(r; R)) = Ael, ref(R) |«'ref(r; R)) , (2) 

which depends parametrically on the nuclear coordinates 
R. Employing an effective single particle description (e.g. 
Hartree-Fock or density functional theory), the electronic 
ground state is given by a single Slater determinant. 


an irrelevant constant, the vibrational part of can be 
rewritten as 

~ 2]\/f ^iic(R) + Eel, ref(R) ~ flaoloQ,, 

a ^ a 

( 5 ) 

where the ladder operators oj,, and aa address the vi¬ 
brational mode OL with frequency The respective di¬ 
mensionless displacement and momentum operators read 

Qa = A ~ Q-a)- 

In the next step, we expand the orbital energies, which 
enter the electronic part of the Hamiltonian, about the 
equilibrium geometry, Rq, of the reference state 
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and similar for the Coulomb interaction 


|vI/,ef(r;R))= n ^-(R)| 0 ). ( 3 ) 

mG{occ.} 

This determinant involves single-particle states (orbitals) 
R)) which are occupied in the reference state. Us¬ 
ing the creation and annihilation operators (R) corre¬ 
sponding to \‘pjn(yi] R)), the electronic part of the Hamil¬ 
tonian can be rewritten 

iJel(R) = Eel, ref(R) + ^ ^m{R){Andm “ ^m) (4) 

m 

+ Y.U^r.{R){dld 

mn 

where em(R) denote the orbital energies. The parame¬ 
ters Sm distinguish between single-particle states that are 
occupied {Sm = 1) or unoccupied {Sm = 0) in the refer¬ 
ence state R)). The energy of a molecular state 

that differs in the population of two electronic orbitals 
m and n from the reference state is not only modified by 
the energy differences em(R) and en(R), but also by the 
Coulomb interaction between the electrons Umn{R)- In 
principle, there is also the necessity for interaction terms 
describing the change in energy for molecular sates that 
differ in the population of three or more electrons from 
the ground state. They are, however, beyond the scope 
of our present considerations. Moreover, in the deriva¬ 
tion of Eq. O, we have applied the adiabatic approxi¬ 
mation, that is neglecting the coupling between different 
electronic states due to electronic-vibrational or electron- 
electron interaction a^^i^^ . Such effects are considered, for 
example, in Refs. [H, HI], For simplicity, we fur¬ 

thermore, do not consider the spin explicitely. 

To characterize the nuclear (vibrational) degrees of 
freedom, we employ the normal modes of the reference 
state |\E'ref(r;R)). The corresponding potential energy 
surface (PES) is given by I4uc(R) + Eei, ref(R). Using 
the harmonic approximation for the PES and dropping 
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Notice that the form invariance of the Coulomb interac¬ 
tion is broken in 0, due to the fact that [/mn{R) rep¬ 
resents an effective Coulomb interaction in a restricted 
model Hamiltonian. As a result, the Hamiltonian used 
to describe the molecular bridge is given by 
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It includes electron-electron interactions Umn , electronic- 
vibrational coupling Xma, and vibrationally dependent 
electron-electron interactions Wmna- 

While the influence of electron-electron interac¬ 
tions and electronic-vibrational coupling on transport 
in molecular junctions have been studied in detail 
beforei*^, vibrationally dependent electron-electron in¬ 
teractions have so far only been considered in different 
contexts^2iS. We first discuss their physical origin. The 
PESs of a generic one dimensional model for a molecu¬ 
lar junction is depicted in Fig. [H where the solid black 
line is associated with the ground state of the neutral 
molecule, the solid red line with the ground state of the 
respective anion, the solid blue line with the first ex¬ 
cited state of the anion and the solid purple line with 
the ground state of the dianion. The minima of these 
PESs represent the equilibrium geometry of the nuclei, 
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FIG. 1. (Color online) Scheme of potential energy surfaces of 
the molecular bridge corresponding to the ground state (GS) 
and first excited state (ES) of the neutral molecule, its anion 
and the dianion. The equilibrium geometry of the respective 
molecular states are marked by dashed vertical lines of the 
same color, the displacement of the equilibrium geometry with 
respect to the ground state of the neutral molecule is given 
by black arrows. 

highlighted by dashed vertical lines. The equilibrium ge¬ 
ometry of the nuclei changes upon charging or excitation 
of the molecule. This constitutes one source of electronic- 
vibrational coupling. For example, the shift of the equi¬ 
librium geometry of the ground state of the anion with 
respect to the ground state of the neutral molecule is 
given by AQi = 2Aii/r2i. Similarly, the equilibrium ge¬ 
ometry of the first excited state of the anion is shifted 
by AQ 2 = 2 A 2 i/ni with respect to the ground state of 
the neutral molecule. Without vibrationally dependent 
electron-electron interactions, the equilibrium position 
of the nuclei in the ground state of the dianion would 
be shifted by AQ 12 = AQi -|- AQ 2 = 2(Aii -|- A 2 i)/ni 
with respect to the ground state of the neutral molecule, 
i.e. would be fixed by the parameters of the PES of the 
singly charged molecule. However, the equilibrium po¬ 
sition of the nuclei of the dianion is not necessarily cor¬ 
rectly described by this shift. The actual nuclear dis¬ 
placement can be characterized by an additional parame¬ 
ter via AQ 12 = AQi-\-AQ 2 + 2Wi2i/^i. This shows that 
vibrationally dependent electron-electron interaction ac¬ 
counts for additional shifts of the equilibrium positions of 
the nuclei that occur if a state differs by the occupation 
of more than one single-particle state from the reference 
state [e.g. |Tref(r;R))). 

The molecule in a molecular junction is coupled to two 
electrodes, which we model by non-interacting electrons, 

■^L/R = X! ^kclck, (8) 

feeL/R 

that are located on the left (L) and the right (R) elec¬ 


trode, respectively. The molecule-lead coupling can be 
described by 

^SL/SR “ ^ ^ + H.C., (9) 

mGM,fcGL/R 

which allows for electron exchange processes between the 
molecule and the leads and determines the level-width 
function 

rL/Ry(e) = 27r ^ ek)- (10) 

feGL/R 

Throughout this article, the leads are modeled as semi¬ 
infinite tight binding chains with inter-site coupling 
strength /3. The corresponding level-width function is 
given by 

rL/Rb(e) = (11) 

with X = e—AtL/R and the Heaviside step function 0Sr^. 
Similar to Vik , the parameters J^L/Ri describe the cou¬ 
pling between the molecular state i and the respective 
lead. For simplicity, we neglect a dependency of the 
molecule-lead coupling on the nuclear degrees of freedom 
or the charge state of the molecule and assume, further¬ 
more, a symmetric drop of bias voltage = —MR = d>/2. 

The overall Hamiltonian used to describe a single¬ 
molecule contact is given by 

H = Ffs + i^L + EfR + Rsl + i?SR. (12) 

B. Effective population-dependent 
electronic-vibrational coupling 

To understand the basic mechanism of vibrationally 
dependent electron-electron interactions, it is expedient 
to diagonalize the molecular part of the Hamiltonian 
Hs- To this end, we employ a generalized small polaron 
transformation^^— of the Hamiltonian 

H = e^He-^ (13) 

with 
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The transformed Hamiltonian H has no explicit 
electronic-vibrational or vibrationally dependent 
electron-electron interaction terms. The effect of these 
interactions is subsumed in the polaron-shifted energy 
levels 
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In a similar way as the bare electronic-vibrational inter¬ 
action Xma influences the electron-electron interaction 
{Umn ^ Umn 2 ^^(Ayja/H(^)), the vibrationally 
dependent electron-electron interaction leads to coupling 
terms proportional to c\cic^^cjc\ck ii ^ j ^ k) and 
clciCjCjclckclci (i yf j yf fc yf /) in H. To be consistent 
with the derivation of Eq. o, we neglect those higher- 
order terms in Eqs. m- 

Aside from an additional renormalization of the bare 
electron-electron interaction strengths, Umn Umn, 
which induces no qualitatively new effects compared 
to the case with Wmna = 0, vibrationally depen¬ 
dent electron-electron interactions affect, in particu¬ 
lar, the structure of the shift operators Xm- Instead 
of c-numbered electronic-vibrational coupling strengths, 
Xma, the shift operators Xm involve effective electronic- 
vibrational coupling strengths 

Xma — Xma 4" ^ ^ ^^mnaid^dji Sn) (^9) 

n^m 

that include the electronic occupation operators {dl^dn — 
Sn) with respect to all single-particle states but the mth 
one. As a result, the effective electronic-vibrational 
coupling strengths Xma depends on the population of 
the single-particle states and thus on the charge of 
the molecule. Specifically, an electron that is trans¬ 
ferred from the electrode to the mth state of the neutral 
molecule, couples with Xma = Xma to the vibrational 
mode a. For an electron that populates state m of the 
charged molecule, where the occupation of the nth state 
differs from the neutral molecule, the effective electron- 
vibrational interaction is Xma = Xma + Wmna- This 
implies the existence of different electronic-vibrational 
coupling strengths for the same electronic state, in the 
above example Xma and Xma + Wmna- The coupling 


strength relevant for an electron entering/leaving the 
molecule depends on the exact population of all the 
other electronic states at the very moment the trans¬ 
port process takes place. An interpretation based on an 
averaged electronic-vibrational coupling strength Xma + 
Wmna {d\dn — 5n) does therefore not provide a correct 
description of the physics. A quantitative description re¬ 
quires an explicit calculation of the transport character¬ 
istics including vibrationally dependent electron-electron 
interactions, which will be the focus of Sec. MM 


C. Master equation approach 

We simulate the transport properties of a single¬ 
molecule junction employing the well established 
Born-Markov master equation methodologyiii^iiSliZSrii. 
Thereby, the central object is the reduced density ma¬ 
trix p, which is obtained as the stationary solution of the 
equation of motion 

dtp{t) = -i[Hs,p{t)\ (20) 

dTtrB{[HsB, [d7sB(T),p(t)pB]]}, 

with 

77 sb ( t ) = e-dHs+Hs)Tjj^^^z(Hs+Hs)T ( 21 ) 

and pB being the equilibrium density matrix of the leads. 
Eq. (1^ can be obtained, e.g. using a second-order ex¬ 
pansion of the exact Nakajima-Zwanzig equatio n^^d^ in 
the coupling 77 sb, including the so-called Markov ap¬ 
proximation. In the applications considered below, we 
will focus on the regime of resonant transport and weak 
molecule-lead coupling, where the Born-Markov master 
equation provides a correct description of the dominating 
transport processes. In this regime, where the molecule 
changes its charge state, the effects of vibrationally de¬ 
pendent electron-electron interaction are also expected 
to be most pronounced. 

We evaluate the master equation (EUl) for the steady 
state, focusing on model systems with two electronic 
states and a single vibrational mode, which is sufficient 
to show the generic effects of vibrationally dependent 
electron-electron interactions. In the calculations, we use 
basis functions \nin 2 )\v) (ni,n 2 S {0,1}) that represent 
the subspace of the electronic |nin 2 ) and the vibrational 
degrees of freedom \i/) {v £ Nq) in occupation number 
representations, respectively. The coefficients of the re¬ 
duced density matrix are thus denoted by 

PZ'n2,n[n', = V>2) 

= {nin2\{iy\p\U)\n[n2). ( 22 ) 

The principal value terms in Eq. (1201) describe the renor¬ 
malization of the molecular energy levels due to the cou¬ 
pling between the bridge and the leads^i^i^l. The im¬ 
portance of these terms has been investigated by Hartle 
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and Millis in a recent study of charge-transfer dynam¬ 
ics in a double quantum dot system^i. For the systems 
considered in this work, where the single particle levels 
are well separated, these terms can be neglected. For the 
same reason we neglect vibrational coherences in the den¬ 
sity matrix, which is a valid assumption for the systems 
considered here, which do not exhibit quasi-degeneracies 
and where the broadening due to molecule-lead coupling 
is small compared to the vibrational energies^i^. 


D. Observables of interest 


For characterizing electron transport through a single¬ 
molecule junction, we analyze the electric current and 
the average vibrational excitation as a function of applied 
bias voltage 

Within the density matrix methodology outlined 
above, the expectation value of an observable O can be 
calculated as 

{0)=tT{pO}= ^ {nin2\{iy\pO\i^) \n 1 n 2 ) (23) 

ni ,n2,i^ 


Specifically, the excitation of the vibrational mode z/ is 
given byM 

{a\,ay) = 

+ - 5j){dldk - Sk))j^ 

jk 

+2 E - SMdk - Sk))jj 

i,j<k ^ 

+2 E ii4d, - d^){d]d, - SMdk - 4) 

i<j,k<l ^ 

x{d\di — 5i))jj (24) 


Thereby, the subscript H/H denotes the Hamiltonian, 
which is used to evaluate the respective expectation 
value. 

The number of electrons entering or leaving the lead 
K {K € {L,R}) per unit time determines the electronic 
current. 


Ik 


{Ik)h = —2e 


dt 


k&K 


(25) 


= 2ie 


_ km 


J2Kk{dlxlck)Tr 


Here, the constant (—e) denotes the electron charge and 
the factor 2 accounts for spin-degeneracy. To second or¬ 
der in the molecule-lead coupling, the current through 
lead K can be written 

poo 

/ dTtrs+B{pSB(T),pBp]-fic}- (26) 
It is noted that this expression is current conserving, i.e. 

Il = —Ir = I- 


III. RESULTS 

To analyze the effect of vibrationally dependent 
electron-electron interactions, we consider a series of min¬ 
imal models comprising two electronics levels and a sin¬ 
gle vibrational mode. The parameters of the models 
are summarized in Tab. HI where, for the sake of clar¬ 
ity, we dropped all vibrational indices and the electronic 
indices of U and W. It is noted that all systems in¬ 
vestigated exhibit relatively weak electronic-vibrational 
coupling strengths. 

The model parameters have been selected to study dif¬ 
ferent aspects and mechanisms, including the manifesta¬ 
tion of vibrationally dependent electron-electron inter¬ 
actions as an effective population-dependent electronic- 
vibrational coupling and its influence on vibrational ex¬ 
citation. Furthermore, signatures of negative differen¬ 
tial resistance and asymmetries in the gate voltage de¬ 
pendences of the current, even in symmetrically coupled 
molecular junctions, are investigated. 


A. Effective population-dependent 
electronic-vibrational coupling 

We start the analysis of the influence of vibrationally 
dependent electron-electron interaction by studying the 
transport properties of the generic model system EFF 
(see Tab.m for different interaction strengths W. In this 
model, the electronic energy levels of the anion and the 
dianion are well separated, which allows for a better iden¬ 
tification of the effect of vibrationally dependent electron- 
electron interaction. As Eq. (ED shows, vibrationally 
dependent electron-electron interaction renormalizes the 
Coulomb interaction strength, U —>■ U, leading to a shift 
of the electronic resonances in the current-voltage char¬ 
acteristics. For better comparison, the value of U is fixed 
and the bare Coulomb interaction strength U is adjusted 
such that the location of the electronic resonances coin¬ 
cides for all values of W. 

Fig. [5] shows the current and the vibrational excitation 
as function of bias voltage $ for model system EFF for 
W = ±0.05 eV. For comparison, also data obtained with¬ 
out vibrationally dependent electron-electron interaction, 
i.e. W = 0 eV, are depicted. The results show that the 
vibrationally dependent electron-electron interaction has 
a significant effect for voltages $ > 2e2. Depending on 
the sign of the interaction W and the specific voltage, it 
results in a decrease or increase of the current compared 
to the system without interaction. Moreover, vibrational 
excitation is enhanced for W = —0.05 eV and reduced 
for W = 0.05 eV. 

These findings can be explained employing the effective 
electronic-vibrational coupling strength Ai introduced in 
Eq. (HID. We start by considering the case W = ±0.05 
eV. Eor low bias voltages, $ ^ 2(ei ±C/), the current and 
the vibrational excitation agree with the noninteracting 
model. At these voltages, the features in the current- 
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TABLE I. Parameters of the model systems investigated in this article. For all calculations, the temperature is set to T = 10 
K. All parameters are given in eV. 


Model 

ei 

^2 

V-L 1 

vb. 1 

VB 2 

vb. 2 

P 

Kl 

Ai 

A 2 

U 

IF 

EFF 

0.15 

0.3 

0.1 

0.1 

0.1 

0.1 

3 

0.1 

0.05 

0.05 

0.525 

0,±0.05 

STMSETUP 

0.15 

0.2 

0.1 

0.01 

0.1 

0.01 

3 

0.1 

-0.05 

0.05 

1.05 

0,±0.05 

DARKST 

0.15 

0.2 

0.1 

0.1 

0.01 

0.01 

3 

0.1 

-0.05 

0.05 

0.05 

0,±0.05 

ASYMM 

0.15 

0.3 

0.1 

0.03 

0.1 

0.03 

3 

0.1 

0.05 

0.05 

0.525 

0,±0.05 




FIG. 2. (Color online) (a) Current-voltage characteristics 
for the model system EFF. (b) Vibrational excitation as a 
function of bias voltage for the system EFF. The vertical 
dashed lines in both plots mark the onset of resonant trans¬ 
port through the corresponding molecular electronic states. 


voltage and vibrational excitation characteristics are re¬ 
lated to the opening of transport channels at energies 
ei/2 + mf2 with m € Z. Populating both electronic states, 
ei and € 2 , is not possible at these low voltages because the 
electrons coming from the leads do not have enough en¬ 
ergy to overcome the Coulomb repulsion. Accordingly, an 


electron impinging on the molecular bridge from the left 
electrode encounters a neutral molecule. Therefore, in 
this elementary charge transport step, the electron cou¬ 
ples to the vibrational degrees of freedom with the cou¬ 
pling strengths Ai = Aid-FF-O = Aj, which is independent 
of the vibrationally dependent electron-electron interac¬ 
tion W. As a consequence, the transport properties a 
low voltages are virtually identical to that of the system 
without the vibrationally dependent electron-electron in¬ 
teraction. 

For higher voltages, $ > 2(ei -b U), new features in 
the current and the vibrational excitation characteristic 
appear which are associated with the opening of trans¬ 
port channels at energies f-1/2 + U + mO with m G Z 
and correspond to transport channels where the electron 
impinging from the left electrode onto the molecule en¬ 
counters a singly occupied molecular bridge. As a con¬ 
sequence, in the elementary charge transport steps, the 
electron couples to the vibrational degrees of freedom 
Ai = Xi+W -1 = \i-\-W. For IF = -1-0.05 eV, the effective 
electronic-vibrational coupling is thus increased, which 
leads to two effects observed in Fig. [2] First, resonant 
transport processes associated with the absorption of vi¬ 
brational energy are enhanced, resulting in an increased 
current for voltages below the onset of resonant trans¬ 
port involving the dianionic molecule, $ < 2 (e 2 + CA).— 
These processes also lead to a decreased vibrational ex¬ 
citation seen in the corresponding voltage regime in Fig. 
Hb). Second, in the weak electronic-vibrational cou¬ 
pling regime considered here, an increased electronic- 
vibrational interaction gives rise to a decreased current 
and a diminished vibrational excitation beyond the onset 
of resonant transport involving the respective electronic 
stateiii^i^^— , in this case the dianionic resonance, at 
voltages $ > 2(e2 + U). 

It is important to note that the current at these higher 
voltages includes both transport channels which couple 
with Ai = Ai to the vibrations and which became ac¬ 
tive already at low bias voltages, and transport chan¬ 
nels that couple with Ai = Ai -b IF. Although the 
average populations of the two electronic levels satu¬ 
rate at about {d[di) « {d\d 2 ) « 0.5 at high bias volt¬ 
ages, we want to stress that the above discussed trans¬ 
port behavior is not described correctly by considering 
transport with electronic-vibrational coupling strengths 
of Ai = Ai -b IF ■ 0.5 corresponding to the average pop¬ 
ulation, but needs to consider the populations of the el- 
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ementary charge transport steps, which are zero or one, 
as discussed above. 

Next we consider the system with W = —0.05 eV. For 
low voltages, as discussed above, the transport is deter¬ 
mined by the effective coupling Xi = Xi + W - 0 and thus 
independent on W. At higher voltages, ^ > 2{ei + U), 
however, the features in the current and the vibrational 
excitation are associated predominantly with transport 
processes, where an electron coming from the left elec¬ 
trode encounters a singly occupied molecular bridge re¬ 
sulting in an effective coupling of Ai = Ai -I- bF ■ 1 = 0 eV 
for the chosen parameters of the model, A = —W = 0.05 
eV. This corresponds to a system where the electronic- 
vibrational coupling and the vibrationally dependent 
electron-electron interaction cancel each other. Com¬ 
pared to the system with W = 0 eV, this results in a 
smaller current for 2ei < e<i> < 2(77 + U) and a larger 
current for e$ > 2(77 + U). At higher bias voltages the 
Franck-Condon blockade of the current is lifted, while 
the smaller current at lower voltages is due to the ab¬ 
sence of transport processes associated with the absorp¬ 
tion of vibrational energy. This is also reflected in the vi¬ 
brational excitation characteristics (Fig. [2Kb)), which for 
W = —0.05 eV exhibits a steady increase with voltage in¬ 
dicating the absence of processes that absorb vibrational 
energy. This increase of vibrational excitation is fur¬ 
ther enhanced because a decreased electronic-vibrational 
coupling leads to an increased vibrational excitation in 
the regime of weak electronic-vibrational coupling due 
to missing electron-hole pair creation processes that ef¬ 
fectively cool the molecular bridge^. As a result, the 
largest average vibrational excitation is observed for the 
system W = —0.05 eV. 


B. Effect of vibrationally dependent 
electron-electron interactions in molecular junctions 
with left-right asymmetry 

Additional effect of vibrationally dependent electron- 
electron interactions arise in models with asymmetric 
coupling to left and right leads as is common, e.g., in 
STM setups. As an example, we consider the model sys¬ 
tem STMSETUP with parameters listed in Tab. HI Due 
to the difference in the coupling to the leads, the elec¬ 
tronic states of the molecular bridge are completely oc¬ 
cupied for positive bias voltages above the onset of res¬ 
onant transport through the respective levels. For neg¬ 
ative bias voltages, they are unoccupied. The anionic 
molecule can be in its electronic ground (7i occupied, 
72 unoccupied) or first excited state (7i unoccupied, 72 
occupied). Depending on the electronic configuration of 
the anion, the dianion is obtained by populating the first 
or the second electronic state. For the specific choice of 
the parameters, Ai = —A 2 = ±kF, one of these channels 
for generating the dianion decouples from the vibrations 
(Ai = 0 eV or A 2 = 0 eV, see below). This allows for 
a better identification of the effect of vibrationally de¬ 
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FIG. 3. (Color online) (a) Current-voltage characteristics for 
the model system STMSETUP. (b) Vibrational excitation as 
a function of bias voltage for the system STMSETUP. The 
vertical dashed lines in both plots mark the onset of res¬ 
onant transport through the corresponding molecular elec¬ 
tronic states. The inset in panel (b) shows an enlargement of 
the region $ « 2(7i/2 -b U). 


pendent electron-electron interaction in the current and 
vibrational excitation characteristics. 

The current-voltage characteristics of model STM¬ 
SETUP, depicted in Fig. [3] (a), exhibits a pronounced 
asymmetry with respect to bias polarity, which is a well 
known effect associated with the asymmetric coupling to 
the leads and the resulting dependence of the electronic 
population on bias polarity^ but barely vary with W. 
However, as shown in Fig. |3Kb), there is a significant in¬ 
fluence of the vibrationally dependent electron-electron 
interaction on the vibrational excitation for positive bias 
voltage. Two features are noteworthy: First, for FF 0 
eV, there is a sudden decrease in vibrational excitation 
at voltages <f> = 2(7i/2 + U), where the dianion becomes 
energetically accessible (see inset of Fig. Elb)). Second, 
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for high bias voltages, the vibrational excitation is signifi¬ 
cantly smaller compared to the W = Q case, independent 
of the sign of W. 

These findings can be explained considering the popu¬ 
lation of the electronic states and the effective electronic- 
vibrational couplings Ai. Due to the asymmetric coupling 
to the leads, the molecule is mostly in the dianionic state, 
once the applied positive bias allows for double charging 
of the molecule. As a consequence, electrons populating 
the molecule from the left encounter an anionic molecule 
most of the times and transport processes with coupling 
^ 1/2 = Ai /2 + bb -1 to the vibrational degrees of free¬ 
dom are dominant. Because one the effective couplings 
vanishes (Ai = 0 eV or A 2 = 0 eV), the average vibra¬ 
tional excitation strongly decreases at the bias voltage 
where transport with A 1/2 = 0 eV becomes energetically 
possible. 

For large positive bias voltages, beyond the onset of 
transport involving the dianion, both molecular elec¬ 
tronic states are almost always occupied. Accordingly, 
the current flowing across the molecule in this bias regime 
is determined by the transport processes that rely on 
the generation the dianion. For the system without vi- 
brationally dependent electron-electron interaction, the 
transport through both electronic states couples with A 
to the nuclear degrees of freedom. Thus both trans¬ 
port channels cause a heating of the vibrational mode. 
For the systems with W ^ 0 eV, on the other hand, 
only the transport channel through one of the electronic 
states couples to the vibrations with twice the coupling 
strength, 2A, whereas the other decouples from the vi¬ 
brations. As a stronger electronic-vibrational coupling 
leads to a decreased vibrational excitation in the regime 
of overall weak electronic-vibrational coupling (i.e. A < 
the systems with W = ±0.05 eV exhibit on av¬ 
erage a smaller number of vibrational quanta for large 
positive bias voltages. 


C. Negative differential resistance 

Another effect introduced by vibrationally dependent 
electron-electron interaction is negative differential resis¬ 
tance (NDR), that is a decrease in current upon increase 
of bias voltage. To demonstrate this effect, we consider 
model system DARKST as specified in Tab. HI It consists 
of a strongly coupled electronic state ei, which is mainly 
responsible for the current flowing through the molecule, 
and a weaker coupled, or dark state € 2 , which mainly 
influences the transport properties of the main channel 
via vibrationally dependent electron-electron interaction. 
Again, the electron-electron interaction has been adapted 
such that the location of the electronic resonances coin¬ 
cides for the systems with IT = 0, ±0.05 eV. 

The current-voltage characteristics for the model sys¬ 
tem, depicted in Fig. SI exhibit one large step at $ = 2ei, 
marking the onset of resonant transport through ei , and 
additional, smaller vibrational features, which depend 
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FIG. 4. (Color online) Current-voltage characteristics of 
model system DARKST. The vertical dashed lines mark 
the onset of resonant transport through molecular electronic 
states. 


on the magnitude and sign of W. The influence of the 
second transport channel £2 is barely visible due to its 
weak coupling to the leads. For the case W = 0, the 
main steps of the current are at $ = 2(ei ± nfl) and 
$ = 2(ei + U + nn) with n € Nq. For W = ±0.05 eV, 
however, there are also distinct features at $ = 2 (e 2 +nfi) 
and $ = 2(e2 + U + nfi). While for W = 0.05 eV the 
current is overall larger than for the reference system 
W = 0, for IT = —0.05 eV it is reduced. Further¬ 
more, for W = —0.05 eV, the vibrational features give 
rise to distinct decreases in the current upon increasing 
bias voltage. This NDR effect is marked by blue arrows 
in Fig. m 

We first consider the results for W = 0.05 eV. In this 
case, the effective couplings for transport through level 
ei are Ai = -0.05 eV and Ai = -0.05 eV +W = 0 eV, 
respectively, depending on the occupation of the levels. 
For transport through level €2 the effective electronic- 
vibrational couplings are A 2 = 0.05 eV and A 2 = 0.05 eV 
+W = O.I eV. Notice that the latter corresponds to an 
enhanced coupling for transport involving the dianion. 
With the population of €2 at biases T > 2e2, the current 
flowing through ei comprises transport channels that in¬ 
volve the dianion. As this transport path decouples from 
the vibrations, Ai = 0 eV, the current is increased com¬ 
pared to transport through ei of the anionic molecule 
with Ai = —0.05 eV. 

Next we study the current-voltage characteristics for 
W = —0.05 eV. In this case, the effective electronic- 
vibrational couplings for level ei are Ai = —0.05 eV 
and Ai = —0.05 eV +W = —0.1 eV. Notice that trans¬ 
port through ei involving the dianion corresponds to an 
enhanced electronic-vibrational coupling. For transport 
through level £ 2 , the effective electronic-vibrational cou¬ 
plings are A 2 = 0.05 eV and A 2 = 0.05 eV +W = 0 eV. 
For $ = 2(ei ± nfi) and d) = 2(ei ± {7 ± nfi) we ob¬ 
serve an increase in the current that is associated with 
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the opening of new transport channels directly populat¬ 
ing ei- For $ = 2(^2 + nfl) and $ = 2(^2 + U + nfl), on 
the other hand, new transport channels directly populat¬ 
ing level 62 become available. Consequently, transport 
processes involving the dianionic molecule become more 
important, thus increasing the significance of transport 
through ei with enhanced electronic-vibrational coupling 
strength Ai = —0.1 eV. As the transport through Ii dom¬ 
inates the current and an enhanced electronic-vibrational 
coupling results in a smaller current, a stepwise decrease 
of the current is observed at bias voltages $ = 2 {e 2 + nQ) 
and = 2 (e 2 -1-17-1- nfl). Notice that this NDR effect 
results from the influence of the level 62 on the trans¬ 
port through ei via the population dependent electronic- 
vibrational coupling Ai. The NDR effect is therefore 
qualitatively different from similar NDR effects such as 
blocking state scenarios^. 

Finally we want to remark that there are model sys¬ 
tems, where the renormalization of the electron-electron 
interaction strength U ^ U caused by the vibrationally 
dependent electron-electron interaction can also result in 
NDR. However, this effect does not appear in the model 
DARKST considered above. 


D. Asymmetries with respect to bias and gate 
voltage 

The vibrationally dependent electron-electron interac¬ 
tion leads to vibrational effects that are influenced by 
the electronic population of the molecular levels. As a 
consequence, vibrational features can change with bias 
polarity and gate voltage. In this section, we study these 
dependences based on conductance maps, that is the con¬ 
ductance as a function of bias and gate voltage. Thereby 
we assume that the only effect of a gate voltage <i>gate 
on the system is to shift the electronic energies of the 
noninteracting molecule, 6 i ^ ei + d)gat(j^. Other inves¬ 
tigations of molecular junction transport based on con¬ 
ductance maps, also referred to as Coulomb diamonds or 
stability diagrams can be found, e.g., 

We start with model system EFF from Sec. IIII A[ 
which is characterized by a symmetric coupling to the 
leads and vibrations. Fig. [5] shows the conductance for 
coupling strengths of IF = 0, ±0.05 eV. The conduc¬ 
tance maps are dominated by the darker red lines, form¬ 
ing a diamond-shaped square in the center. These lines 
correspond to the onset of resonant transport through 
the electronic levels. Since the systems are corrected 
for the energy shift introduced by the vibrationally de¬ 
pendent electron-electron interaction, these features ap¬ 
pear at the same positions for any value of IF. Addi¬ 
tionally, the conductance maps exhibit a distinct struc¬ 
ture of less pronounced lines, which are associated with 
the onset of vibrationally coupled transport and differ 
in magnitude with IF. For IF = —0.05 eA^, fewer vi¬ 
brational lines are visible and for IF = ±0.05 eV, the 
vibrational structures are observed over a wider range 


of voltages compared to the IF = 0 eV case. For 
IF = 0 eV, the conductance maps exhibits three symme¬ 
tries with respect to the transformations $bias —> —‘hbias, 
‘hgate —!■ —‘hgate —(ei±e2±17) and consequently also with 
respect to $bias, $gate -^bias, “^gate “ (d ± £2 ± U). 
The first is due to the symmetric coupling of the molecule 
to the leads and is unaffected by the vibrationally depen¬ 
dent electron-electron interaction. The second and third 
symmetry are broken for IF 7 ^ 0 eV. Lastly, the data 
show blue patterns, corresponding to small negative dif¬ 
ferential resistance, located at high bias voltages but also 
in areas close to the onset of resonant transport. 

These findings can be rationalized as follows. The loca¬ 
tion of the lines related to electronic and vibronic trans¬ 
port are unchanged by IF because the vibrationally de¬ 
pendent electron-electron interaction does not alter the 
vibrational energy. For IF = —0.05 eV, fewer vibrational 
transport channels exist as transport involving the dian¬ 
ionic resonance effectively decouples from the nuclear de¬ 
grees of freedom. For IF = ±0.05 eV, transport including 
the dianionic molecule couples with twice the strength to 
the nuclear displacement such that processes including 
several vibrational quanta are more pronounced than in 
the IF = 0 eV case. The symmetry with respect to the 
gate voltage is lifted because the vibrational effects de¬ 
pend on the electronic population of the molecule. This 
leads to a change of the effective electronic-vibrational 
coupling upon variation in gate voltage, and hence to 
a different vibrational structure. The origin of the NDR 
for high bias voltages is the finite bandwidth of the leads, 
which are modeled as semi-infinite chains. The NDR for 
the model system with IF = 0.00 eV for gate voltages 
around —0.2 V, —0.475 V and —0.7 V is caused by the 
influence of the vibrational nonequilibrium state on the 
transport properties of the junction as discussed in Ref. 
[ 2 ^. The change in the NDR structure for IF 7 ^ 0 eV is re¬ 
lated to the change in the effective electronic-vibrational 
coupling for transport involving the dianion, resulting in 
an altered vibrational nonequilibrium state and the effect 
discussed in Sec. imm 

Next, we consider the model system ASYMM with 
asymmetric molecule-lead coupling, which allows to 
study the relation between the vibrationally dependent 
electron-electron interaction and bias polarity. The cor¬ 
responding conductance maps for IF = 0, ±0.05 eV are 
shown in Fig. [5] As a consequence of its asymmetric cou¬ 
pling to the leads, the electronic population of this model 
is sensitive to the polarity of the applied bias and so 
is the influence of the vibrationally dependent electron- 
electron interaction. As in the previous model, we ob¬ 
serve a change in the vibrational structure of the con¬ 
ductance map upon a variation in IF. For IF = —0.05 
eV, the lines corresponding to vibrational transport are 
less pronounced, in particular for gate voltages below 
Ri —0.45 V only few vibrational features are present. 
For IF = ±0.05 eV, the vibrational features are more 
pronounced and are observed for all gate voltages. The 
data also exhibit weak NDR effects, which are influenced 
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FIG. 5. (Color online) Conductance as a function of bias and gate voltage for the system EFF for W — 0, ±0.05 eV. Notice 
the different scale used for negative conductances. 
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FIG. 6. (Color online) Conductance as a function of bias and gate voltage for the system ASYMM for W = 0, ±0.05 eV. Notice 
the different scaling for negative conductances. 


by the vibrationally dependent electron-electron interac¬ 
tion. Due to the asymmetric coupling to the leads, there 
is no symmetry in the conductance maps with respect 
to bias polarity, $bias — ^bias- Remarkably, also the 
symmetry of the conductance map under the transfor¬ 
mation ^biasj ^gate ^ ^biasi ^gate (ci ± 62 ± ff), is 
broken for W ^ 0. In the VF = 0 eV case, this symme¬ 
try is a consequence of the fact that both bias and gate 
voltage change the average population of the electronic 
states in a similar way and that there is no distinction 
between transport involving the anion or the dianion. In¬ 
cluding vibrationally dependent electron-electron inter¬ 
action, transport involving the anion or the dianion is no 
longer equivalent. As a result the transport mechanism 
depends on the total charge of the molecule and, there¬ 
fore, this symmetry of the conductance map is broken. It 
is interesting to note that most of the conductance maps 
measured in experiment display asymmetrie a^d^i^^~ i^ 


IV. CONCLUSION 

We have investigated the effect of vibrationally de¬ 
pendent electron-electron interactions in single-molecule 
junctions. The additional interaction is a result of the 
dependence of the Coulomb interaction on the nuclear 
displacement and accounts for the fact that vibronic 
transport processes depend on the charge state of the 
molecule, i.e. are different for transport through, e.g., an 
anionic or dianionic state of the molecule. Employing a 
generalized small polaron transform, we have shown that 
vibrationally dependent electron-electron interaction re¬ 
sults in an effective electronic-vibrational coupling, which 
depends on the electronic occupation, and can be used 
to rationalize the effects on charge transport. 

Employing a master equation approach, we have an¬ 
alyzed the basic mechanism and the manifestations of 
vibrationally dependent electron-electron interactions in 
single-molecule junctions. Depending on the strength an 
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the sign of the interaction it may result in a significant 
alteration of the transport characteristics and the vibra¬ 
tional nonequilibrium excitation and may cause NDR. 
For selected values of the interaction strength, the in¬ 
terplay between electronic-vibrational interaction and vi- 
brationally dependent electron-electron interaction can 
also lead to regimes where electronic-vibrational coupling 
is effectively switched off. In junctions with asymmetric 
molecule-lead coupling, vibrationally dependent electron- 
electron interaction may cause a strong dependence of 
the vibrational nonequilibrium excitation on the bias po¬ 
larity. Finally, vibrationally dependent electron-electron 
interaction can give rise to asymmetries of conductance 
maps related to the different description of transport in¬ 
volving the anionic and the dianionic molecule. The lat¬ 
ter finding may be of particular interest in the context of 


experimental results, which have the tendency to display 
asymmetrie s^' ^ 
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